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I The critical phenomenon of the zero temperature superfluid— Bose-glass phase 

, transition for hard-core bosons on a three-dimensional disordered lattice is studied 

Q I using a quantum real-space renormalization-group method. The correlation-length 
O ' 

exponent i' and the dynamic exponent z are computed. The critical exponent z 

^ I is found to be 2.5 for compressible states and 1.3 for incompressible states. The 



exponent is shown to be insensitive to z as that in the two-dimensional case, and 
has value roughly equal to 1. 
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Disordered boson system provides a prototypical example of zero temperature quantum 
critical phenomena Experimentally, such a disordered boson system can be realized in 
liquid He^ in random media , or in disordered superconductors |Q where cooper pairs can 
be modeled as composite bosons. As the amount of disorder is varied, these systems exhibit 
a continuous phase transition from the superfiuid (SF) phase to a disordered (Bose-glass 
(BG)) phase. Understanding of this quantum critical phenomenon has been a subject of a 
considerable amount of recent experimental and theoretical studies. 

Besides a diverging length scale (the correlation length) ^, the SF-BG transition is also 
characterized by a diverging time scale r. Denoting 6 as the distance to the criticality, ^ and 
T can be described by the critical exponents u and z, defined by ^ oc 6^" and r oc oc 6~'^^. 
From a general scaling argument. Fisher et al. concluded that for compressible states z 
is equal to d, the dimensionality of the system, while z = 1 for systems with long-range 
Coulumb interactions in any dimension 0. The critical exponent u cannot be deducted 
directly from the scaling theory, but a rigorous lower bound has been established, i.e., 
z/ > 2/(i 1^. Unfortunately, standard field-theoretical renormalization-group (RG) method, 
which is proven to be a powerful approach in the study of critical phenomena, has been eluded 
so far from being applied to this SF-BG phase transition. A major difficulty associated with 
this approach is due to the lack of a proper zero loop (mean field) theory describing the SF- 
BG transition at finite dimension 0, upon which perturbation series, such as e-expansions, 
can be developed. For this reason, a different RG formulation, namely the real-space RG 
(RSRG) approach, becomes a useful alternative for investigating the critical phenomena in 
disordered boson systems. 

RSRG has been applied previously to one- dimensional (ID) and two-dimensional 

(2D) systems ||^. In this paper, we apply the RSRG method developed previously by 
Zhang and Ma ^ to investigate the SF-BG transition in three dimension. Such a study is 
important not only because of its intrinsic theoretical interest, but also due to its relevance 
to experiments. Three dimensional (3D) disordered boson systems are directly realized 
in liquid He^ in Vycor, aerogels, xerogels, and other random media [^. While present 
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experiments are mainly focused on the effect of disorder on tlie finite temperature superfluid 
pfiase transition, refined experiments are hopefully to be capable of extracting information 
about the zero temperature critical point in the near future. It is clearly desirable to perform 
calculations for critical exponents in 3D which characterize the physical properties of the 
system in the vicinity of the transition. 

In the following we shall briefly outline the RSRG scheme which has been presented in 
detail in Ref. 0. The RSRG scheme yields, when applied to the ID systems, no (non- 
trivial) fixed point, indicating the instability of the superfluid phase against any amount 
of disorder for hard-core boson systems, in agreement with other ID RG calculations 
and exact results |10|. For 2D and 3D systems, it gives a non-trivial fixed point separating 
the superfluid phase and the disordered phase. In 2D, the critical exponent z was found 
to be about 1.7 for compressible states and about 0.9 for incompressible states. Since z 
is close to 2 and 1 in these two cases, and since it increases slowly with increasing block 
size 0], this calculation can be viewed as a supporting evidence of the scaling prediction 
by Fisher et al. in 2D. The critical exponent v is found to be quite insensitive to the 
type of states, and is roughly equal to 1.4 [[^. This value of v satisfies the rigorous lower 
bound of Chayes et al. H, and is also consistent with other numerical estimates O]. The 



product zv for incompressible states may be directly compared with data obtained from 
experiments on superconductor-insulator phase transitions in homogeneous amorphous films 
0J^. For 3D systems, on the other hand, no systematic computation has yet been performed 
to our knowledge. Other numerical methods used for ID and 2D systems, such as exact 



diagonalization [|TT| or quantum Monte Carlo simulations |jT2|] , seem quite formidable for 3D 
disordered systems within the present computational capacity. Thus, the RSRG method 
provides a unique way to probe the physical properties of the 3D SF-BG transition. 

The system under consideration is a lattice model of hard-core bosons with random 
potential [|1|, 

n = -tY^iblb, + H.c.) + Y^iw, - ^)h]h, , (1) 



where fo] and bj denote boson creation and annihilation operators at lattice site j. 
indicates nearest-neighbor summation and Wj is the random on-site potential with (inde- 
pendent) Gaussian distribution. The hard-core constraint is enforced by the requirement 
that at each site the occupation number b^jbj equals to either or 1. While it is clearly a 
simplification to the realistic systems, this model is believed to have captured the essential 
physics of the zero-temperature SF-BG phase transition [jl],0. 

This hard-core boson model is equivalent to a quantum spin-| XY-mode\ with transverse 
random fields 

n = -JY: is^sj + sfs]) - y: h,s; , (2) 

<i,j> j 

via the mapping Sj ^ 6j, Sj' ^ bj, J ^ 2t, and hj ^ ^, — Wj. These two equivalent 
representations of the problem provide a convenient way to study the underlying physics of 
the system and we will use them alternatively throughout our discussions. 



Our real-space RG method consists of the following steps: 

i) break the lattice into blocks of size n^; 

ii) compute the block spin which is given by two low energy eigenstates of the block 
Hamiltonian; Since = J2j Sj (corresponding to the particle number Np = J2j b]bj in 
the boson language) is a good quantum number, eigenstates of H are also simultaneously 
eigenstates of (Np). As described in Ref. 0, this can be accomplished either by (a) 
selecting the lowest states of two chosen subspaces of particle number q and g + 1, or by (6) 
selecting the two lowest states among the ground states of the block Hamiltonian for each 
subspace of definite number of particles. The two states chosen in this way are found to have 
adjacent particle number q' and g' + 1. The field acting on the block spin is given by the 
energy difference between the two states. In general, the block fields will follow a different 
distribution than that of site fields. We chose to keep track only the mean h = hi and the 
variance h = {hf — hi'Y^'^ of the renormalized field, and thus map it onto a Gaussian. 

Hi) calculate the effective couplings between the block spins, which is given by the 
nearest-neighbor couplings of the site spins between two adjacent blocks; Due to the pres- 



ence of disorder, couplings between block spins are also randomized in the RG iterations. 
However, these block couplings are found to be always positive, and the renormalized system 
remains unfrustrated. This allows one to approximate the block couplings by its mean, and 
thus confine the RG iterations within the original parameter space. 

iv) repeat the RG iteration defined above to find fixed point (s) and compute critical 
exponents. 

Here we compute the exponents z and u for 3D systems. For computational details, 
see Ref. [0. In the procedure (a), fluctuations in particle number at any given region is 
one. Since the block size is fixed and it allows only a discrete set number of particles 
occupying a block, states described by such a procedure are incompressible. The procedure 
(6), which is valid only at the 'particle- hole symmetrical' filling (half boson per site in our 
case) where density is conserved through statistical fluctuations in the random fields [0, 
yields a compressible ground state. 

In the present work we used a simple-cubic lattice with block size n<j = 2 x 2 x 2. 
Gaussian random fields were generated numerically and typically an average of 3 x 10^ 
random configurations was performed. In the previous work [0, calculations for 3D systems 
were carried out only for the incompressible state with q = 6. Here we complete the study 
for the incompressible states with other q values, and investigate the compressible states 
as well. Results for fixed points and critical exponents are summarized in the table. For 
incompressible states with other q values, they can be obtained through the 'particle- hole 
symmetry' 0. 

As anticipated, critical exponents obtained for incompressible states at different density 
have roughly the same value, indicating that the same universality class is probed. As that in 
the 2D case, the critical exponent u in 3D is shown to be rather insensitive to the procedures 
adopted, roughly equal to 1.0. The critical exponent z, on the other hand, is quite different 
for the compressible and incompressible states. The present RSRG calculation yields z = 2.5 
and 1.3 for the compressible and incompressible states, respectively. 

Our results for u in both compressible and incompressible states satisfy the lower bound 
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of Chayes et al. , and provide the first systematic calculation on critical exponents for 3D 
systems. General scaling argument suggests that z = d ii the system is compressible 0. 
For systems with long range (Coulomb) force, dimensional analysis indicates that z = 1 m 
any spatial dimension . If one assumes that they belong to the same universality class as 
those incompressible states considered here, this value should be compared with our results 
obtained through the fixed q procedure. Our RSRG results for z in 3D deviate more from 
these scaling predictions than those obtained for the 2D case 0. The value of z deviates 
from the scaling result by 17% for the compressible states, and by 30% for the incompressible 
cases. The corresponding deviation for 2D systems is 15 and 10%, respectively. It is worth 
to remark that within our RG procedure, the value of z for compressible states differs from 
that for incompressible states roughly by a factor of two in both 2D and 3D cases. While 
in 2D this is exactly what the scaling argument suggested, it is clearly not the case in 3D. 
Our computation is limited to the cubic block of size = 8, thus we cannot address the 
question of size effects. The next isotropic block has the size = 3'^ = 27 for which the 
calculation is computationally not feasible for the time being. Thus we are not able to test 
the direction in which the result would converge with increasing block size. Since the RSRG 
calculation is expected to approach the exact result as the block size increases |^, such a 
finite-size study is highly desirable. On the other hand, from the study of 2D systems one 
expects that the RSRG estimate for critical exponents itself is rather insensitive to small 
changes of the block size. 

Although the RG procedure (b) for the compressible states is valid only at the 'particle- 
hole symmetric' point where particle number is conserved through the statistical fluctuation 
in random fields, it nevertheless allows one to probe the parameter space away from this 
special point by studying the flow of the RG iterations. The RG flow diagram is similar to 
what one gets in 2D with even block size Ug 0. Again, the (unstable) flxed point at /i = 
axis controls the critical phenomena of the SF-BG phase transition. 

Before closing, we briefly comment about the RSRG approach to disordered quantum 
systems. Since the block states are computed through an isolated block Hamiltonian, one 
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would worry about the effect of the (long-range) coherent quantum fluctuations which may 
not be properly taken into account. A procedure to overcome this difficulty by incorporating 
different boundary conditions for the block states has been recently suggested, and has been 
applied to noninteracting fermion systems [0. While a direct application of such a RG 
method to disordered system has not been possible within present computational facilities, 
we would like to remark that by averaging over random configurations one can partly achieve 
the goal that one wishes to accomplish by averaging over the boundary conditions. Indeed, it 
seems quite ironic that our RSRG procedure works precisely in the 'strong' disordered regime 
where the effect of fluctuations can be incorporated through averaging over the randomness, 
and losses its validity as the pure limit is approached 0. 

In summary, we have studied a hard-core disordered boson system in 3D cubic lattice 
using a quantum RSRG method. We have found that there exists a non-trivial fixed point 
describing the zero temperature SF-BG phase transition. The critical exponent z for in- 
compressible states is about 1.3, and about 2.5 for compressible states. The exponent v 
is insensitive to and is roughly equal to 1.0. As in the 2D case, only one universality 
class is found, and the critical behavior of the zero temperature SF-BG phase transition is 
controlled by the 'particle-hole symmetric' fixed point. 
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TABLES 



state {h/J)* ^ -2 

incompressible [q = Q)\ 1.6 1.0 1.2 

incompressible (g = 5) 2.4 1.0 1.3 

incompressible (g = 4) 2.8 1.0 1.3 

compressible 6.0 1.1 2.5 



t calculated in ref. 0]. 
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